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^ ■ Abstract 

(^ ■ We consider QCD radiative corrections to Standard Model Higgs boson production 

in association with a W boson in hadron collisions. We present a fully exclusive 
calculation up to next-to-next-to-leading order (NNLO) in QCD perturbation theory. 
To perform this NNLO computation, we use a recently proposed version of the 

/S ■ subtraction formalism. Our calculation includes finite-width effects, the leptonic decay 

c^ ' of the W boson with its spin correlations, and the decay of the Higgs boson into a bb 

pair. We present selected numerical results at the Tevatron and the LHC. 
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The search for the Higgs boson or its alternatives in theories beyond the Standard Model (SM) 
is the key endeavour of current high-energy colliders. The Tevatron has recently excluded at the 
95% CL a SM Higgs boson in the mass range 158 GeV < m^ < 173 GeV [Ij. At the LHC the 
Higgs search has now started. If the LHC will confirm the Tevatron exclusion, the attention will 
move to the low mass region, where the Higgs search is more difficult. 

One of the most important production mechanisms of a light Higgs boson at hadron colliders 
is the so called Higgs-strahlung process, i.e. the associated production of the Higgs boson together 
with a vector boson V {V = W^, Z). At the Tevatron, this is the main search channel in the 
low mass region, mn ^ 140 GeV, since the lepton(s) from the decay of the vector boson provide 
the necessary background rejection. At the LHC, this production channel has been considered 
less promising, due to the large backgrounds. Recent studies [2J have indicated that at large 
transverse momenta, employing modern jet reconstruction and decomposition techniques, WH 
and ZH production can be recovered as promising search modes for a light Higgs boson at the 
LHC. 

In order to fully exploit this channel it is important to have accurate theoretical predictions 
for the production cross section and the associated distributions. Theoretical predictions with 
high precision demand in turn detailed computations of radiative corrections. The next-to-leading 
order (NLO) QCD corrections to V77 production are the same as those of the Drell-Yan process [3]. 
Beyond NLO, the QCD corrections differ from those to the Drell-Yan process [1] by contributions 
where the Higgs boson couples to the gluons through a heavy-quark looptO. The impact of these 
additional terms is, however, expected to be rather small. Using the classical result of Ref. [1], 
the computation of the NNLO inclusive cross section has been completed in Ref. 0. Also NLO 
electro-weak corrections have been evaluated [6j. 

The effect of NNLO QCD radiative corrections on the inclusive cross section is relatively 
modest [5]. However, it is important to study how QCD corrections impact the accepted cross 
section and the relevant kinematical distributions. This is particularly true when severe selection 
cuts are applied, as it typically happens in Higgs boson searches. 

The evaluation of higher-order QCD radiative corrections to hard-scattering processes is well 
known to be a hard task. The presence of infrared singularities at intermediate stages of the 
calculation does not allow a straightforward implementation of numerical techniques. In particular, 
fully differential calculations at the NNLO involve a substantial amount of conceptual and technical 
complications [7]-[l2]. In e^e~ collisions, NNLO differential cross sections are known only for 
two [131 [E] and three jet production [121111]. In hadron-hadron collisions, fully differential cross 
sections have been computed only in the cases of Higgs production by gluon fusion [TTJ [TSl [12] 
and of the Drell-Yan process [201I2H I22]. 

In this Letter we focus on WH production and we present the fully exclusive NNLO QCD 
computation for this process. The calculation is performed by using the subtraction formalism of 
Ref. [18] , and it is based on an extension of the numerical program of Ref. [21] . We include finite- 
width effects, the leptonic decay of the W boson, with its spin correlations, and the decay of the 
Higgs boson into a bb pair. Only diagrams in which the Higgs boson is radiated by a VT boson are 
considered, i.e. we neglect the contributions in which the Higgs boson couples to a heavy-quark 



'I'ln the case of ZH production there are also gluon-gluon induced terms where both the Z and H bosons couple 
to a heavy-quark loop [5] . 



loop. Comparing our results to those of Ref. |23], where NLO predictions for WH + jet, including 
these additional diagrams, are presented, we expect the neglected contributions to be at the 1% 
level or smaller. When no cuts are applied, and the W and H are produced on shell, our numerical 
results agree with those obtained with the program VHONNLO [5]. 

In the following we present an illustrative selection of numerical results for WH production 
at the Tevatron and the LHC. We consider u, d, s, c, b quarks in the initial state and we use 
the (unitarity constrained) CKM matrix elements Kd = 0.97428, K. = 0.2253, Kb = 0.00347, 
Vcd = 0.2252, Vcs = 0.97345, Kb = 0.0410 from the PDG 2010 [24J. As for the electroweak 
couplings, we use the so called G^ scheme, where the input parameters are Gp, niz, mw In 
particular we use the values Gp = 1.16637 x 10"^ GeV-^, mz = 91.1876 GeV, mw = 80.399 GeV 
and Tiv = 2.085 GeV. Throughout the paper, we consider a SM Higgs boson with mass itih = 120 
GeV and width Th = 3.47 MeV [2S]. We compute the H ^ bb decay at tree level in the massless 
approximation, and we normalize the Hbb Yukawa coupling such that BR{H — )■ bb) = 0.649 |25j . 
We use the MSTW2008 [26] sets of parton distributions, with densities and as evaluated at each 
corresponding order (i.e., we use (n + l)-loop as at N"LO, with n = 0, 1,2). The central values 
of the renormalization and factorization scales are fixed to the value /i/j = fip = tt^w + ^^h- 

We start the presentation of our results by considering WH production at the Tevatron {pp 
collisions at a/s = 1.96 TeV). We use the following cuts (see e.g. Ref. [27j). We require the charged 
lepton to have transverse momentum p^ > 20 GeV and pseudorapidity \rii\ < 2, and the missing 
transverse momentum of the event to fulfil p^^^^ > 20 GeV. Jets are reconstructed with the kp 
algorithm with R = 0.4 [28]. We require exactly two jets with pp > 20 GeV and \ri\ < 2, and at 
least one of them has to be a 6 jet, with |?7| < 1. 



a(fb) 


LO 


NLO 


NNLO 


/iF = /i/j = (mw + mH)/2 


4.266 ±0.003 


4.840 ± 0.005 


4.788 ±0.013 


llF = l^R = mw + TTlH 


3.930 ±0.003 


4.808 ± 0.004 


4.871 ±0.013 


fJ'F = fJ-R = 2{mw + mn) 


3.639 ±0.002 


4.738 ± 0.004 


4.908 ±0.010 



Table 1: Cross sections for pp — > WH ± X — )■ luibb + X at the Tevatron. The applied cuts are 
described in the text. 

In Tabled] we report the accepted cross sectioro at LO, NLO and NNLO in the case of three 
different scale choices around the central value fip = ^J'R = f^w + ^h- The impact of NLO 
(NNLO) corrections ranges from +13% to ±30% (—1% to ±4%), depending on the scale choice. 
The scale dependence is at the level of about ±1% both at NLO and NNLO. 

In Fig. [1] we show the transverse-momentum spectrum of the dijet system at LO, NLO, NNLO. 
The lower panel of the figure shows the ratio NNLO/NLO. We see that the shape of the spectrum 
is rather stable, when going from NLO to NNLO, within the statistical uncertainties. 

We now consider VFi7 production at the LHC (y^ = 14 TeV). We follow the selection strategy 
of Ref. [2] (see also [29]): the Higgs boson is searched for at large transverse momenta through its 
decay into a collimated bb pair. We require the charged lepton to have pip > 30 GeV and \rii\ < 2.5, 



■l-Throughout the paper, the errors on the vahies of the cross sections and the error bars in the plots refer to an 
estimate of the numerical errors in the Monte Carlo integration. 



1.50 

1.25 

1.00 

0.75 

0.50 

0.25 

0.00 
1.1 
1.0 
0.9 
0.8 



- pp ^WH+X^li/bb +X 



^ 



Vs=1.96 TeV 

m„=130 GeV 
MF=MR=mw+mH 

o-Lo = 3.930 ± 0.003 fb 
CTni,o = 4.808 ± 0.004 fb 
<^NNL0=4.871 ± 0.013 fb 




100 



200 
Pt" (GeV) 



300 



400 



Figure 1: Transverse-momentum spectrum of the dijet system for pp — )■ WH + X — )■ lubb + X at 
the Tevatron at LO (dots), NLO (dashes) and NNLO (solid). The applied cuts are described in 
the text. 

and the missing transverse momentum of the event to fulfil p™'^^ > 30 GeV. We also require the 
W boson to have p^ > 200 GeV. Jets are reconstructed with the Cambridge/Aachen algorithm 
[30] . with R = 1.2. One of the jets (fat jet) must have p^ > 200 GeV and \rij\ < 2.5 and must 
contain the bb pair. There should not be other jets with pt > 20 GeV and \ri\ < 5. 



a(fb) 


LO 


NLO 


NNLO 


fiF = I^R = {mw + mH)/2 


2.640 ± 0.002 


1.275 ±0.003 


1.193 ±0.017 


I^F = I^R = mw + ruH 


2.617 ±0.003 


1.487 ±0.003 


1.263 ±0.014 


jJ'F = I^R = 2{mw + mH) 


2.584 ±0.003 


1.663 ±0.002 


1.346 ±0.013 



Table 2: Cross sections for pp — )■ WH + X — )■ luibb + X at the LHC. The applied cuts are described 
in the text. 

In Table [2] we report the corresponding accepted cross sections at LO, NLO and NNLO. We 
see that the impact of NLO and NNLO corrections is negative and larger than at the Tevatron. 
The NLO (NNLO) effect ranges from -52% to -36% (-6% to -19%), depending on the scale 
choice. The scale dependence goes from about ±13% at NLO to about ±6% at NNLO. The NLO 
and NNLO accepted cross sections are compatible within the scale uncertainties. 

In Fig. [2] we show the transverse-momentum spectrum of the fat jet at LO, NLO and NNLO. 
The lower panel shows the ratio NNLO/NLO. We see that the shape of the distribution is relatively 
stable when going from NLO to NNLO. 



We add few comments on the different impact of QCD radiative corrections at the Tevatron 
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Figure 2: Transverse-momentum spectra of the fat jet for pp — )■ WH + X — )■ lubb + X at the LHC 
at LO (dots), NLO (dashes) and NNLO (solid). The applied cuts are described in the text. 



and at the LHC [31J. At the Tevatron, the invariant mass of the VKi/ system is M^h ~ mw + '^'H- 
The typical scale of the accompanying QCD radiation is of the order of about (1 — z) Mi^u where 
(1 — 2;) = (1 — M'^^/s) is the average distance from the partonic threshold. The effect of the veto 
on additional jets is thus marginal if the jet veto scale, p™*°, is of the order of (1 — z) Mwh- In this 
case the perturbative expansion appears under good control. The situation at the LHC is different 
in two respects. First, the invariant mass of the VFi/ system is larger, due to the high px required 
for the W and the Higgs candidate. Second, the typical distance from the partonic threshold is 
larger, i.e. {1 — z) is larger than at the Tevatron, due to the increased -y/i. As a consequence, 
a stringent veto on the radiation recoiling against the WH system spoils the cancellation of the 
infrared singularities in real and virtual corrections, and contributions enhanced by the logarithm 
of the ratio (1 — z)Mwh/p^^° are definitely relevant. We have checked that the reduction in the 
accepted cross section is in fact due to the jet veto, the impact of QCD corrections being positive 
if the jet veto is removed. 

We have illustrated a calculation of the NNLO cross section for WH production in hadron 
collisions. The calculation is implemented in a parton level event generator and allows us to apply 
arbitrary kinematical cuts on the W and H decay products as well as on the accompanying QCD 
radiation. We have studied the impact of NNLO QCD corrections in two typical cases at the 
Tevatron and the LHC. At the Tevatron, the perturbative expansion appears under good control. 
At the LHC, by searching for events where the Higgs boson is boosted at high p^, the impact of 
QCD corrections is more sizeable, and the stability of the fixed-order calculation is challenged. 
More detailed studies, along the lines of Refs. [321 [33] , are needed in order to assess the relevance 
of these fixed-order perturbative result. 
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